Cosmological Signatures of the Interaction 
between Dark-Energy and Massive Neutrinos 
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We investigate whether interaction between massive neutrinos and quintessence scalar field is the 
origin of the late time accelerated expansion of the universe. We present cosmological perturbation 
theory in neutrinos probe interacting dark-energy models, and calculate cosmic microwave back- 
ground anisotropies and matter power spectrum. In these models, the evolution of the mass of 
neutrinos is determined by the quintessence scalar field, which is responsible for the cosmic acceler- 
ation today. We consider several types of scalar field potentials and put constraints on the coupling 
parameter between neutrinos and dark energy. Assuming the flatness of the universe, the constraint 
we can derive from the current observation is m " < 0.87eV at the 95 % confidence level for the 
sum over three species of neutrinos. We also discuss on the stability issue of the our model and on 
the impact of the scattering term in Boltzmann equation from the mass- varying neutrinos. 

PACS numbers: 98.80.-k,98.80.Jk,98.80.Cq 



Interoduction: After SNIaQ and WMAP0 obser- 
vations during last decade, the discovery of the acceler- 
ated expansion of the universe is a major challenge of 
particle physics and cosmology. In order to understand 
unknown 76% components of the critical density of the 
universe with a negative pressure (dark-energy) , the pos- 
itive cosmological constant term seems to be a serious 
candidate for the dark energy. In this case the cosmo- 
logical constant A and it's energy density A/8irG remain 
constant with time and correspnding mass density p\ — 
6.44 x 10- 30 (Q A /0.7)(/i/0.7) gem- 3 , where h is the Hub- 
ble constant Hq expressed in units of 100 kms~ 1 Mpc~ 1 
and Q\ = 0.76. Although cold dark matter model with a 
positive cosmological constant Q {AC DM) provides an 
excellent explanation of the SNla data, the present value 
of A is 10 123 times smaller than the value predicted by 
the particle physics model. Among many alternative can- 
didates for dark energy^, H, H| the scalar field model like 
quintessence is a simple model with time dependent w, 
which is generally larger than — 1. Because the differ- 
ent w leads to a different expansion history of the uni- 
verse, the geometrical measurements of cosmic expansion 
through observations of SNla, CMB, and Baryon Acous- 
tic Oscillations (BAO) can give us tight constraints on 
w. Further, if the dark energy is dynamical component 
like a scalar field, it should carry its density fluctuations. 
Thus, the probes of density fluctuations near the present 
epoch, such as cross correlation studies of the integrated 
Sachs- Wolfe effect 0, H| and the power of Large Scale 
Structure (LSS) @, can also provide useful information 
to discriminate between cosmological constant and oth- 
ers. Yet, current observational data can give only poor 
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constraints on the properties of dark energy fluctuations 
[ifl [Tl| . Another interesting way to study the scalar field 
dark energy models is to investigate the coupling between 
the dark energy and the other matter fields. In fact, a 
number of models which realize the interaction between 
dark energy and dark matter, or even visible matters, 
have been proposed so far [lj|, 0, 0, dM [l6[- Obser- 
vations of the effects of these interactions will offer an 
unique opportunity to detect a cosmological scalar field 

Interacting Dark-Energy with Neutrinos: In 

this letter we investigate the cosmological implication 
of an idea of the dark-energy interacting with neutrinos 
[Til [li| . For simplicity, we consider the case that dark- 
energy and neutrinos are coupled such that the mass of 
the neutrinos is a function of the scalar field which drives 
the late time accelerated expansion of the universe. In 
previous works [H, [l!| , potential term was treated as a 
dynamical cosmology constant, which can be applicable 
for the dynamics near present epoch, but kinetic energy 
term has been ignored in thier discussions. However the 
kinetic contributions become important to describe cos- 
mological perturbations in early stage of universe, which 
is fully considered in our work. 

Equations for quintessence scalar field are given by 
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where V{4>) is the potential of quintessence scalar field, 
Vi(4>) is additional potential due to the coupling to neu- 
trino particles [l9l . [20l ] , and m„ (</>) is the mass of neutrino 
couDled to the scalar field. TL is hi a. where the dot reD- 
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resents the derivative with respect to the conformal time 

T. 

Here we consider three different types of the 
quintessence potential: (1) inverse power law potentials 
(Model I), (2) SUGRA type potential models (Model II), 
(3) exponential type potentials (Model III), which are 
given, respectively: 

M 4 ( MEL)* M 4 ( MsiYe 3 *'^ : mV^'. 



where 
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Energy densities of mass varying neutrino (MVN) and 
quintessence scalar field arc described as 
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From equations ^ and ([7]l. the equation of motion for 
the background energy density of neutrinos is given by 

fi In in 

P„ + 3H(P» + P.) = - 3P„) ■ (9) 



The evolution of neutrinos req uires solving the Boltz- 
mann equations in the case: [2ll |22|: 



dq 1 • 2 m dm dx 1 

Tr =—h l3 qn n>-a - — — 



(10) 



Our analytic formula in eq. (|10[) is different from those of 
[32| and [33| , since they have omitted the contribution of 
the varying neutrino mass term. The first order Boltz- 
mann equations written in the synchronous gauge reads 
HI: 



= -i q -(n.k)kS^ dl ^ fo .(11) 
e q z dip dmq 

The Bolzmann hierarchy for neutrinos, obtained expand- 
ing the perturbation ^ in a Legendre series can be written 
as fUlii : 
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Constrains on the MaVaNs parameters: As was 

discussed in the introduction, the coupling between cos- 
mological neutrinos and dark energy quintessence could 
modify the CMB and matter power spectra significantly. 
It is therefore possible and also important to put con- 
straints on coupling parameters from current observa- 
tions. For this purpose, we use the WMAP3 0, HH] 
and 2dFGRS H§ data sets. 

The flux power spectrum of the Lyman-a forest can 
be used to measure the matter power spectrum at small 
scales around z < 3 [13, It has been shown, however, 
that the resultant constraint on neutrino mass can vary 
significantly from ^m v < 0.2eV to 0.4eV depending on 
the specific Lyman-a analysis used [2i|. The complica- 
tion arises because the result suffers from the systematic 
uncertainty regarding to the model for the intergalactic 
physical effects, i.e., damping wings, ionizing radiation 
fluctuations, galactic winds, and so on [3(j. Therefore, 
we conservatively omit the Lyman-a forest data from our 
analysis. 

Because there are many other cosmological parame- 
ters than the MaVaNu parameters, we follow the Markov 
Chain Monte Carlo(MCMC) global fit approach [3l[ to 
explore the likelihood space and marginalize over the nui- 
sance parameters to obtain the constraint on parame- 
ter (s) we are interested in. Our parameter space consists 
of 



P = (n b h 2 ,Q c h 2 ,H,T,A s 



n,, m,, a 



,/3) , (17) 



where uibh 2 and £l c h 2 are the baryon and CDM densities 
in units of critical density, H is the hubble parameter, r is 
the optical depth of Compton scattering to the last scat- 
tering surface, A s and n s are the amplitude and spectral 
index of primordial density fluctuations, and (mj, a, 0) 
are the parameters of MaVaNs. 



TABLE I: Global analysis data within la deviation for differ- 
ent types of the quintessence potential. 



Quantites 


Model I 


Model II 


Model III 


WMAP-3 data 


Q,b h 2 [W] 
SIcdm h 2 [10 2 ] 


2.21 ±0.07 


2.22 ±0.07 


2.21 ±0.07 


2.23 ±0.07 


11.10 ±0.62 


11.10 ±0.65 


11.10 ±0.63 


12.8 ±0.8 


H 


65.97 ±3.61 


65.37 ±3.41 


65.61 ±3.26 


72 ±8 


Z re 


10.87 ±2.58 


10.89 ± 2.62 


11.07 ± 2.44 




a 


< 2.63 


< 7.78 


< 0.92 






< 0.46 


< 0.47 


< 0.58 




n s 


0.95 ±0.02 


0.95 ±0.02 


0.95 ±0.02 


0.958 ±0.016 


A 3 [W 10 ] 
ft Q [10 2 ] 


20.66 ± 1.31 


20.69 ± 1.32 


20.72 ± 1.24 




68.54 ±4.81 


67.90 ± 4.47 


68.22 ±4.17 


71.6 ±5.5 


Age/Gyrs 

n MV N h 2 [w 2 ] 


13.95 ±0.20 


13.97 ±0.19 


13.69 ±0.19 


13.73 ±0.16 


< 0.44 


< 0.48 


< 0.48 


< 1.97(95%C*.L.) 


r 


0.08 ±0.03 


0.08 ±0.03 


0.09 ±0.03 


0.089 ± 0.030 



Larger (3 will generally lead larger m v in the early uni- 
verse. This means that the effect of neutrinos on the den- 



sity fluctuation of matter becomes larger leading to the 
larger damping of the power at small scales. A complica- 
tion arise because the mass of neutrinos at the transition 
from the ultra-relativistic regime to the non-relativistic 
one is not a monotonic function of (3. Even so, the cou- 
pled neutrinos give larger decrement of small scale power, 
and therefore one can limit the coupling parameter from 
the large scale structure data. 

As shown in table HI we find no observational signa- 
ture which favors the coupling between MaVaNs and 
quintessence scalar field, but we don't need to finetune 
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FIG. 1: Differences between the CMB power spectra with 
and without the scattering term in the geodesic equation of 
neutrinos with the same cosmological parameters. 
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the coupling parameters, and obtain the upper limit on 
the coupling parameter j3 as 

P < 1.11, 1.36, 1.53 (2a), (18) 

and the present mass of neutrinos is also limited to 

fi„/it d ay < 0.0095, 0.0090, 0.0084 (2a), (19) 

for models I, II and III, respectively. 
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FIG. 2: (Left panel): Typical evolution of the effective mass 
of the quintessence scalar field relative to theHubble scale, for 
all models considered in this paper. (Right panel): Typical 
evolution of the sound speed of neutrinos c s = 8P v /8p v with 
the wavenumber k — 2.3 x 10 -3 Mpc -1 , for models as indi- 
cated. The values stay positive stating from 1/3 (relativistic) 
and neutrinos are stable against the density fluctuation. 



Results and discussions: Here we discuss two im- 
portant points of this work and results of our analysis 
: (a) the impact of the scattering term of the Boltzmann 
Equation, (b) the instability issue in our models, and (c) 
neutrino mass bounds in the interacting neutrino dark- 
energy models. 

(a) Impact of the new scattering term: Re- 
cently, perturbation equations for the MaVaNs models 
were nicely presented by Brookfield et al. [32j , (see also 
[33l ]) which are necessary to compute CMB and LSS spec- 
tra. A main difference here from their works is that we 
correctly take into account the scattering term in the 
geodesic equation of neutrinos, which was omitted there 
(see, however, [HI). Because the term is proportional to 
^ and first order quantity in pert urbation, our results 
and those of earlier works [32I |33| remain the same in 
the background evolutions. However, as will be shown 
in the appendix, neglecting this term violates the energy 
momentum conservation law at linear level leading to the 
anomalously large ISW effect. Because the term becomes 
important when neutrinos become massive, the late time 
ISW is mainly affected through the interaction between 
dark energy and neutrinos. Consequently, the differences 
show up at large angular scales. In Fig. (Q]), the differ- 
ences are shown with and without the scattering term. 
The earlv ISW can also be affected bv this term to some 



extent in some massive neutrino models and the height 
of the first acoustic peak could be changed. However, the 
position of the peaks stays almost unchanged because the 
background expansion histories are the same. 

(b) Instability issue: As shown in [H, |3(|, some 
class of models with mass varying neutrinos suffers from 
the adiabatic instability at the first order perturbation 
level. This is caused by an additional force on neutri- 
nos mediated by the quintessence scalar field and oc- 
curs when its effective mass is much larger than the hub- 
ble horizon scale, where the effective mass is defined by 

To remedy this situation one should con- 



m 



cfT 



sider an appropriate quintessential potential which has a 
mass comparable the horizon scale at present, and the 
models considered in this paper are the case [32]. Inter- 
estingly, some authors have found that one can construct 
viable MaVaNs models by choo sing certain couplings 
and/or quintessential potentials 37, [H, [3§|. Some of 
these models even realises m e s 3> H. In Fig.©, masses 
of the scalar field relative to the horizon scale m e s/H 
are plotted. We find that m c g < H for almost all period 
and the models are stable. We also dipict in Fig.@ the 
sound speed of neutrinos defined by c 2 s = 8P v /5p v with 
a wavenumber k — 2.3 x 10~ 3 Mpc -1 . 

(c) Neutrino Mass Bounds: When we apply the re- 
lation between the total sum of the neutrino masses M v 
and their contributions to the encrav densitv of the uni- 
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verse: f2„/i 2 = M„/(93.14eV), we obtain the constraint 
on the total neutrino mass: M v < 0.87eV(95%C.L.) in 
the neutrino probe dark-energy model. The total neu- 
trino mass contributions in the power spectrum is shown 
in Figj3j where we can see the significant deviation from 
observation data in the case of large neutrino masses. 



FIG. 3: Examples of the total mass contributions in the mat- 
ter power spectrum in Model I (Left panel) and Model III 
(Right panel). For both panels we plot the best fitting lines 
(green dashed), lines with larger neutrino masses M v =0.3 
eV (blue dotted) and M v = 1.0 eV (cyan dot-dashed) with 
the other parameters fixed to the best fitting values. Note 
that while lines with M v = 0.3 eV can fit to the data well 
by arranging the other cosmological parameters, lines with 
M v = 1.0 eV can not. 
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